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In  this  paper,  we  describe  the  theoretical  and  experimental  determination 
of  the  integrated  time  dependent  concentration  distribution  of  an 
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We  have  derived  this  result  for  the  case  of  constant  applied  flux 
(chronopotentiometry) .  The  experimental  technique  used  for  observing  the 
integrated  concentration  distribution  over  the  electrode  was  imaging  by  a  low 
noise  charge  coupled  device  camera  with  time  and  spatial  resolution  of 
approximately  20ms  and  2Q@m  respectively.  The  experiment  demonstrates  the 
contribution  of  the  tertiary  current  distribution  across  the  surface  of  the 
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Introduction 


There  has  been  considerable  interest  recently  in  the  behavior  of 
microelectrodes  because  of  their  unique  applications  in  electrochemistry 
(1-15).  Due  to  the  ease  of  fabrication  and  convenience  of  use,  disk  and  ring 
electrodes  are  the  most  commonly  used  geometries.  In  addition  to  measurements 
in  media  of  high  resistance,  small  disk  and  ring  microelectrodes  have 
advantages  for  kinetic  measurements  since  quasi  spherical  diffusion  fields  are 
established  more  rapidly  than  at  line  or  band  electrodes.  These  fields  result 
in  high  rates  of  mass  transport  to  the  surface.  As  one  application, 
therefore,  depending  on  the  size  of  the  electrode,  the  kinetics  of 
heterogeneous  electron  transfer  reactions  and  of  fast  reactions  in  solution 
may  be  studied  by  steady  state  measurements . (16) .  Large  disk  electrodes  (ca. 
1mm  or  larger)  generally  will  not  demonstrate  steady  state  conditions  before 
convection  dominates  the  mass  transport. 

Mathematical  analysis  and  experimental  results  have  Indicated  that  the 
diffusion  limited  flux  for  fast  heterogeneous  reactions  approaches  very  large 
values  at  the  electrode- insulator  boundary  of  embedded  electrodes.  The 
combined  effects  of  finite  rates  of  heterogeneous  reactions  and  non-uniform 
distribution  of  potential  and  current  across  the  electrode  surface  will  limit 
the  total  reaction  rate  at  the  edges  for  real  systems(l). 

While  there  has  been  a  variety  of  analytical  and  simulation  methods 
developed  for  determining  the  complex  transient  electrochemical  responses  at 
finite  disk  electrodes  (2-5,16-29),  only  recently  have  exact  solutions  been 
developed  for  a  variety  of  steady  state  and  transient  techniques.  Recently  - 
reported  the  analyses  of  the  responses  for  a  variety  of  transient 
electrochemical  experiments  based  on  the  properties  of  discontinuous  Integra . 
(1-5,16).  In  this  paper,  we  develop  the  theory  for  the  time  dependent 
integrated  concentration  distribution  of  an  electrogenerated  species  above  . 
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finite  disk  electrode  under  constant  flux  conditions  (chronopotentiometry) 
which  is  valid  at  all  times  for  any  size  circular  disk  electrode.  The 
predicted  time  dependent  distribution  is  compared  to  the  experimental  response 
obtained  from  the  in  situ  imaging  the  surface  of  a  disk  electrode  with  a 
charge  coupled  device  (CCD)  camera  located  above  the  plane  of  the  electrode. 
The  experiment  is  similar  to  that  performed  by  Engstrom  (30),  in  which  an 
image  of  surface  flux  was  obtained  by  photography  at  various  times.  We  extend 
his  results  by  performing  the  time  and  space  dependent  experiments  in  situ, 
and  by  comparing  the  results  to  an  exact  analysis  of  the  constant  flux  model 
of  the  finite  disk  electrode.  It  is  observed  that  the  effects  of  non-uniform 
distribution  of  current  (tertiary  distribution)  across  the  electrode  distorts 
the  predicted  concentration  distribution. 

Experimental 

Concentration  detection  was  made  by  an  optical  camera  system  baaed  upon  a 
charge  coupled  device  (CCD) .  The  detection  system  consisted  of  a  front 
illuminated,  uncoated  Thomson  CSF  CCD  chip  (Model  TH7882CDA,  pixel  size  23/xm  X 
23jim,  physical  size  8.8mm  x  13.2mm)  housed  in  a  cryogenic  camera  head  (Model 
CH210,  Photometries,  Ltd.,  Tucson,  AZ) ,  and  operated  at  -120  C.  Data  were 
acquired  via  the  associated  readout  electronics  (Model  CC200  and  Model  CE200. 
Photometries,  Ltd.,  Tucson,  AZ) ,  operated  under  GPIB  control  by  a  Compaq 
386/20  host  computer.  Timing  of  optical  detection  was  provided  by  an  external 
circuit  driving  the  camera  head  mounted  shutter,  which  was  synchronized  with 
the  electrochemistry.  Simultaneous  recording  of  the  electrochemical  current 
transients  was  accomplished  with  a  Lecroy  9400  Digital  Oscilloscope , and 
transferred  to  the  host  computer  via  the  GPIB.  The  optical  cell  was  a 
standard  UV-VIS  specular  reflectance  configuration,  outfitted  with  a  7mm 
diameter  platinum  working  mirror  electrode  which  was  illuminated  by  a  quartz 
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halogen  lamp.  Wavelength  control  was  accomplished  with  appropriate 
Interference  filters.  The  solution  was  lOmM  benzophenone  In  lOOmM 
tetra-n-butylammonium  tetrafluoroborate  In  dry  acetonitrile. 

Chemicals . 

Acetonitrile  (Aldrich  Gold  Label)  was  used  as  received.  Benzophenone 
(J.T.  Baker  Chemical  Co.,  Phillipsburg,  NJ)  was  recrystallized  from  hexane  and 
dried  under  vacuum.  Tetra-n-butylammonium  tetrafluoroborate  (TBAF, 
Chem-Biochem  Research,  Salt  Lake  City,  UT)  was  used  as  received.  All 
solutions  were  prepared  and  transferred  under  dry  nitrogen. 

Theoretical  Considerations 

General  Time  and  Space  Dependent  Solution 

We  have  previously  solved  the  time  dependent  diffusion  equation  in 
circular  cylindrical  coordinates  for  a  simple  electrochemical  experiment 
involving  a  single  reactant  and  product  (2-5): 

A  +  e'  £  B  [1] 

The  equation  is 


dC  _  D  3^C  +  D  3C  +  d^C 
dt  dr2  r  dr  dzZ 


[2] 


where  C  is  the  concentration  of  the  reactant  or  product,  and  r  is  the  radial, 
distance  coordinate  measured  from  the  center  of  the  disk  electrode.  We  assume 
that  the  disk  electrode  is  imbedded  in  the  insulating  plane  z  -  0.  The 
initial  condition  for  the  experiment  is 
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r  >  0,  z  >  0,  t  -  0,  C  -  C®  [3] 

(where  C®  is  the  bulk  concentration) ,  and  the  boundary  conditions  are 

0  S  r  <  .,  z  -  0.  t  >  0  ,  d[|§]  -  •  Q. 

r  >  a,  z  -  0,  t  >  0.  d[|§]  -  0.  [4] 


for  a  constant  flux  (chronopotentiometric)  experiment.  In  the  previous  work, 
(3),  we  obtained  the  general  result  (in  Laplace  space) 


cp  go 

5  "  i - Ds“  JQ  ®xp[-f(A.q)z]Jo(mr)Ji(aa)f^T  [5] 

where 


[6] 


Here,  s  is  the  Laplace  transformation  variable,  while 


a2-[f(A.q)]2-  q2. 


[7] 


and  J  and  J  are  Bessel  functions  of  the  zeroth  and  first  order  respective l\ 
oi  r 

We  rewrite  [5]  as 


C"i - D~  J  exp[-(q2+  Q2)z]Jo(or)Ji(oa)— [8] 

o  s  (q  +  a  ) 


In  the  experiments  described  in  this  work,  we  record  the  integrated 
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concentration  distribution  over  the  entire  region  above  the  electrode  surface 


0  <  z  < 


We  therefore  integrate  [8]  over  all  z 


flO  cp  ao 

J  C  dz  -  C t  -  -jp-  J  J  exp[-(q2+  a2)z]  J^arJJ^aa) 


dz  da 

,  2,  2. 1/2 
s  (q  +  a  ) 


-  f  t  . 

o  s  (q  +  a  ) 


“  Q*  J  JQ(ar)Ji(aa)- 


s  (s  +•  Da  ) 


We  invert  under  the  integral  sign 


t  - ±— - -  -  1  -  exp 

I  s(s  +  Da2)  J  Da2 


to  obtain 


M} 


^  A  +  ft 


1-exp  -Dazt  da 


1 1-exp  \\dp 
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where 
0  —  at 
and 


[14] 


[15] 


Equation  [13]  is  the  integrated  time  dependent  concentration  distribution 
profile  over  the  disk. 

At  sufficiently  long  times  and  small  values  of  the  flux,  we  do  not 
observe  a  transition  time  and  always  reach  the  steady  state  value  for  a  disk 


(16) 


C 

Av 


r*  8Qa 
C  '  3^D  ' 


[16] 


If  Q  is  sufficiently  large  we  observe  a  transition  time  as  the  surface 
concentration  of  the  reactant  approaches  zero.  For  this  condition  and  for  any 
disk  of  radius  a  we  have  (3) 


from  which  the  transition  time  can  be  obtained.  In  the  experiments  describ.  : 
herein,  we  have  maintained  the  flux  at  sufficiently  low  values  so  as  not  to 
reach  a  transition. 


6 


Results  tad  Discussion 

The  development  of  the  time  dependent  integrated  concentration 
distribution  over  a  circular  disk  electrode  of  electrogenerated  benzophenone 
ketyl  anion  in  a  constant  flux  (chronopotentiometric)  experiment  was  recorded 
as  described  in  the  experimental  section.  Even  though  the  total  flux  to  the 
electrode  is  constant  in  time,  the  flux  is  not  constant  across  the  surface  of 
the  electrode.  At  a  disk  electrode,  a  non-linear  current  distribution  will 
develop  across  the  surface  due  to  the  discontinuity  of  the  electrode  surface 
at  the  electrode- insulator  interface.  At  short  tines,  the  higher  electric 
field  strength  at  the  edge  of  the  disk  compared  to  the  interior  planar  portion 
will  cause  preferential  migration  of  ions  to  the  edge.  The  double  layer  at 
the  edge  is  therefore  charged  at  a  faster  rate  than  at  the  interior  points, 
and  electron  transfer  begins  first  at  the  edge.  Migration  of  supporting 
electrolyte  ions  toward  the  edge  will  continue  due  to  electrogeneration  of 
ionic  species  of  opposite  charge  there.  As  a  result,  there  is  further 
depletion  of  ions  from  regions  over  the  interior  of  the  the  disk  toward  the 
edge  causing  further  delay  in  the  charging  of  the  surface.  This  situation 
leads  initially  to  higher  concentrations  of  ketyl  near  the  edge  of  the 
electrode.  Figure  1  is  a  plot  of  the  data  obtained  for  the  integrated 
absorbance  profiles  at  40,  100,  200,  300,  400,  500,  and  1000ms,  respectively 
The  integrated  distributions  in  the  Figure  are  a  convolution  of  the 
concentration  distributions  predicted  by  equation  [13]  and  those  generated  b. 
the  effects  of  the  tertiary  current  distribution.  Figure  lb  are  the  predict..  . 
responses  according  to  equation  [13]  using  literature  values  for  the  diffusi..:. 
coefficient.  Figure  2  repesents  differences  in  absorbances  for  various  times 
After  approximately  100ms  (determined  by  successive  incremental  subtractions 
we  find  that  the  flux  becomes  essentially  constant  across  the  surface  of  the 
electrode.  When  this  condition  is  attained,  we  predict  and  observe  that  the 
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intensity  of  the  absorbance  of  the  concentration  profiles  increases  linearly 
with  time  across  the  surface,  except  near  the  edge  where  they  to  fall 
asymptotically  toward  the  bulk  value  (Cg-  0)*  At  long  times,  we  predict  that 
the  concentration  distribution  will  asymptotically  decrease  from  the  center  of 
the  disk  since  there  is  no  flux  at  z  -  0,  r  >  a.  This  "flat"  distribution 
will  be  observed  at  shorter  times  for  small  disks  than  at  larger  ones,  Figure 
3,  and  is  a  function  of  the  parameter  a/(Dt)  ,  Equation  [13].  We  have 
pointed  out  that  these  profiles  are  similar  to  those  in  constant  potential 
experiments  in  which  there  is  heterogeneous  kinetic  control  (1).  As  an 
illustration,  the  curve  a/(Dt)  -  0.01  corresponds  to  the  integrated 
concentration  profile  developed  in  a  chronopotentiometric  response  over  a  disk 
microelectrode  of  radius  0.1/im,  1  second  after  application  of  the  constant 
flux  Q  (D  -  1.0  x  10  8cm2s*1),  whereas  the  bottom  curve  (a/(Dt)1/2-10.0) 
corresponds  to  the  integrated  profile  which  has  developed  under  the  same 
experimental  conditions  after  a  microsecond.  Alternately,  for  a  1cm  radius 
electrode,  the  bottom  curve  corresponds  to  a  time  of  10*s  !  For  the 
experiments  performed  herein,  therefore,  we  expect  to  observe  a  linear 
distribution  of  concentration  above  the  electrode  except  very  near  the  edge 
(after  relaxation  of  the  tertiary  current  distribution) . 

The  nature  of  the  experiment  precludes  accurate  determination  of 
integrated  concentration  distributions  outside  r  -  a.  The  reflection  of  light 
is  specular  from  the  electrode  surface,  whereas  reflection  from  the  insulator 
surface  at  r  >  a  is  both  diffuse  and  inherently  less  intense.  An  additional 
problem  can  arise  due  to  the  non-uniform  lapping  of  edges  at  the  metal - 
insulator  boundary.  Inhomogenieties  across  the  boundary  will  always  be 
present  due  to  mechanical  polishing.  These  points  are  much  more  obvious  in 
experiments  involving  surface  flux  (30)  (surface  luminescence)  imaging,  and 
will  be  the  subject  of  another  report  (31).  As  a  result,  the  error  in 
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determining  the  transmission  becomes  limited  by  the  digitization  accuracy  of 
the  instrumentation.  We  have  therefore  truncated  the  data  at  points  r  >  a 
since  they  are  not  significant.  The  observation  of  time  dependent 
concentration  and  surface  flux  distribution  profiles  above  the  surface  of 
small  finite  disk  electrodes  requires  more  sophisticated  focussed  imaging 
techniques  and  multichannel  analysis  techniques.  In  this  context,  it  is 
interesting  to  compare  the  differences  in  the  profiles  for  the  cases  of 
constant  surface  flux  (chronopotentiometry)  and  constant  average  surface 
concentration  (chronoamperometry) .  This  point  and  applications  to  small 
electrode  imaging  will  be  reported  shortly. 
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Figure  legends 


Figure  1 

Figure  2 

Figure  3 


.  (a).  Integrated  concentration  distributions  recorded  (bottom  to  top) 
at  40,  100,  200,  300,  400,  300,  and  1000ms,  respectively,  r  and  a 
are  the  radial  coordinate  and  the  radius  (7mm)  of  the  electrode 
respectively.  The  profiles  are  for  electrogenerated  benzophenone 
ketyl  anion  at  a  platinum  electrode,  (b).  The  predicted  responses 
for  the  same  times  recorded  in  (a) . 

Differences  in  radial  integrated  concentration  distributions  across 
the  disk  electrode  for  the  system  described  in  Figure  1  at  (bottom 
to  top)  400- 300ms,  500- 300ms,  and  1000- 500ms  respectively. 

.  Plots  of  the  dimensionless  integrated  concentration  distribution 
across  the  surface  of  a  disk  electrode  /  insulator  assembly  as  a 

function  of  the  dimensionless  parameter  a/(Dt)  (Equation  [13]). 
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